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We report a combined valence band photoemission and Auger spectroscopy study of single crys- 
talline Ca(Fe,Co)2As2 and Ba(Fe,TM)2As2 with TM=Ni or Cu. The valence band photoemission 
data show directly that the TM-states move to higher binding energies with increasing atomic num- 
ber, contributing less and less to the states close to the Fermi level. Furthermore, the 3d 8 final state 
of the LVV Auger decay, which is observed for Ni and Cu, unambiguously reveals the accumulation 
of charge at these impurities. We also show that the onsite Coulomb interaction on the impurity 
strongly increases when moving from Co over Ni to Cu. Our results quantify the impurity potentials 
and imply that the superconducting state is robust against impurity scattering. 



PACS numbers: 74.70.Xa, 71.20.-b, 79.60.-i, 82.80.Pv 

High-temperature superconductivity (HTS) in the iron 
pnictides is one of the most intensively studied topics in 
current condensed matter science. Although research ef- 
forts established many properties of these materials, im- 
portant questions still remain to be clarified pQ . One of 
these questions regards the effect of replacing a few per- 
cent of Fe by other transition metals (TMs) like Co, Ni or 
Cu. The substitution of Fe by Co or Ni induces HTS [I] 
and therefore it is very important to understand the ef- 
fects of the TM-impurities on the electronic structure. 
Exactly these effects, however, remain controversial (e.g. 
Refs. l3)4l and references therein). 

It is often assumed that Co, Ni and Cu mainly dope 
electrons into the Fe bands without seriously affecting the 
electronic structure otherwise [5] . But this interpretation 
has been challenged recently. While some studies con- 
firm that the TM substitutions indeed dope additional 
charge carriers into the Fe-bands [SHE] , others report the 
lack of additional charge at the Fe sites [9HTT] , casting 
doubts on the electron doping via TM substitution. In 
addition, recent theoretical studies show that the random 
scattering potential introduced by the substituted TMs 
has important effects, which are often neglected but may 
in fact cause significant modifications to the electronic 
system [21 [3j [T2l [13] . Hence, there is a considerable ongo- 
ing debate about the effects of TM substitutions in the 
iron pnictides. 

In this letter we shed light on this issue and investigate 
the effects of Co, Ni and Cu-substitution by means of 
photoemission spectroscopy (PES), Auger electron spec- 
troscopy (AES) as well as model calculations. Specifically 
we report studies of Ca(Fe,Co)2As2 and Ba(Fe,TM)2As2 
(TM=Ni, Cu). For Co we chose the Ca- instead of Ba- 
material, in order to avoid complications due to the over- 



lap of the Ba M^g and the Co £2,3 edges [4|. The exper- 
imental data were collected at the UE52-PGM beamline 
of the HZB-BESSY II synchrotron source. Single crys- 
talline samples were grown as described elsewhere [HI ITS] 
and cleaved in situ just before the measurements, which 
resulted in fresh and mirror-like surfaces. More details 
about the sample characterization can be found in the 
supplementary information. PES and AES spectra were 
recorded using a Scienta R4000 hemispherical analyzer, 
whereas the X-ray absorption (XAS) was measured via 
the total electron and fluorescence yield. The photon en- 
ergy was calibrated against the Au 4/ lines, which means 
that the measured binding energies always refer to the 
chemical potential. All spectra where recorded with the 
photon beam at an incident angle of 35° with respect to 
the sample surface, which was parallel to the FeAs-layers. 
For the PES and AES measurements horizontal and ver- 
tical polarized light was used, respectively, to optimize 
the relative intensities of the two components [161. 

The valence band photoemission data for 
Ca(Fe,Co) 2 As 2 and Ba(Fe,TM) 2 As 2 (TM=Ni, Cu) 
are presented in Fig.[TJ including the spectra of the cor- 
responding pure parent compounds. In the case of the 
parent materials, the PES intensity at 0-2 eV and 3-6 eV 
is due to mostly Fe:3d and As:4p derived bands [T71 118j . 
The comparison of Ca(Fei_ a; Co a ;)2As2 with x — 0.06 
to its parent compound CaFe2As2, reveals that both 
spectra are identical within the error of the experiment. 
However, for the Ni- and Cu-substituted materials, 
additional features can be observed clearly in the valence 
band photoemission. An additional structure centered 
at the binding energy (e^ 4 ) = 1.4 eV is present in the 
Ni-substituted BaFe2As2, which does not exist in the 
parent material and which is evident in the difference of 
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FIG. 1. Valence band photoemission data for 

BaFe 2 As 2 , CaFe 2 As 2 , Ca(Fei_ ;c Co :r )2As2 with x = 0.06, 
Ba(Fei_ I Ni^)2As2 with x = 0.05, and Ba(Fei_ a; Cu a ;)2As2 
with x = 0.09. All measurements were performed at room 
temperature using a photon energy of hp — 600 eV and 
normalized by the area above of 7 eV binding energy. 

the two spectra. Even more striking, a strong additional 
peak at (eg") = 4eV is found for the Cu-substituted 
material. We stress that these additional features were 
fully reproducible and observed for different samples. 

The PES data in Fig{T] is experimental proof that the 
3d impurity states move to higher binding energy with 
increasing atomic number Z. While the Ni states around 
1.4 eV are still located inside the Fe bands, although close 
to the band bottom, the Cu impurity states are well be- 
low the Fe bands. The lack of additional features for the 
Co-substituted material can be attributed to the fact that 
these states hybridize strongly with the Fe bands and are 
distributed over essentially the same energy region. 

The above results are in good agreement with pre- 
vious reports, where similar trends have also been ex- 
tracted |4| [13] . An important new observation here is 
that replacing Fe by heavier Ni and Cu, removes Fe- 
states close to the Fermi level and adds impurity states 
at higher binding energy, as is nicely demonstrated by 
the data in Fig.[l] The PES data thereby shows that 



the heavier TMs contribute less to the states close to the 
Fermi level. Note that these data can only provide quali- 
tative information about changes in the density of states 
(DOS), because the PES cross sections for Fe, Co, Ni and 
Cu are different and not known precisely. 

Resonant PES data for the three substituted samples 
are presented in Fig. [2] including the XAS and the non- 
resonant PES. The intensity maps as a function of photon 
energy (hu) and binding energy (eg) have been obtained 
by subtracting the non-resonant PES spectrum, in order 
to highlight the resonant features. In all cases a strong 
increase of the PES intensity can be observed, as soon 
as the photon energy is tuned to the TM-impurity L3 
edge. The PES maps in Fig. [2] show the typical Auger 
features of resonant PES, namely (i) a resonant Raman 
Auger decay right at the TM L 3 threshold, which oc- 
curs at a fixed binding energy, and (ii) the conventional 
Auger decay at higher hu, which leads to a constant ki- 
netic energy of the emitted Auger electron (e^ s ) and 
therefore shifts linearly when plotted as a function of 
e B = hv- e£ZS Comparing the PES maps for 

Co, Ni and Cu, it can immediately be observed that the 
Auger final states move to higher binding energies with 
increasing Z and, at the same time, develop a clear fine 
structure. 

We will consider the following processes: 3d" + hu — > 
2p3d n + e* -> 3d' 1 - 2 + k + e* , with 2p the 2p core hole, 
e* the corresponding photo-excited electron, and k the 
emitted Auger electron. The first step corresponds to the 
absorption of the incoming photon and the second step is 
the LVV Auger decay. Note that e* does not participate 
in the LVV Auger decay, resulting in 2 additional holes 
in the 3d valence shell. In the case of an isolated atom, 
the interaction between these holes splits the Auger final 
state into a characteristic multiplet, consisting of differ- 
ent terms with different energy, total spin S and angular 
momentum L. Each of these final LS'-terms gives rise 
to a line in the total spectrum with a specific intensity. 
If these lines are sharp and well defined, the multiplet 
structure of the AES provides a unique fingerprint of the 
electronic final state configuration. 

The broad AES spectrum of Co, however, does not 
show any clear multiplet structure, as can be observed 
in Figs.[2]a) and d). This implies that the Co-states are 
strongly hybridized with the Fe-host: the state with the 
2 holes on Co couples strongly to the continuum of states 
where one hole has moved from Co into the Fe-host, re- 
sulting in the observed broad and structureless spectrum. 
This observation is in agreement with a recent resonant 
PES study on Ca(Fei_ :E Co :E )2As2 [2"3] . 

As can be seen in Figs.[2]b) and e), Ni shows a much 
more structured LVV spectrum. The observed lineshape 
is essentially identical to that found for Ni- metal [16] [20] , 
showing structures that agree very well with an atomic 
Ni:3<i 8 Auger final state multiplet. Both the energy split- 
ting between the different terms and the relative inten- 
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FIG. 2. Resonant PES scans across the L 3 edge for a) Ca(Fe,Co)2As2, b) Ba(Fe,Ni) 2 As 2 , and c) Ba(Fe,Cu)2As 2 . On the left 
side of each map in the top row the XAS spectra are displayed. From each map a off-resonance PES spectra was subtracted, 
which is shown below each map. The AES peak for each incident photon energy is indicated by black dots, showing the 
crossover from a resonant Raman (e_g =const.) to the normal AES decay (ck =const.) . The lower panels show LVV Auger 
spectra corresponding to the cuts indicated by the dashed lines in the upper panels. Term splittings of the atomic 3d 8 final 
state were calculated using Cowan's code [TU] and the AES intensities were taken from Refs. 1201 and 1211 (vertical lines). The 
calculations do not take into account the full resonance process, but consider the Auger decay of the core hole state only. 



sities obtained by the atomic calculation (20] are in very- 
good agreement with the experimental data. 

The observation of a 3d 8 Auger final state for Ni 
is remarkable, because it is incompatible with the of- 
ten assumed electron doping. The latter would result 
in a Ni:3d 6 in the ground state, from which the 3c? 8 
Auger final state cannot be reached at all. A Ni:3e? 6 
is hence excluded by our data. Instead, the data im- 
plies a strong Ni:3d 9 component in the ground state. 
Here the 3d 8 final state can be reached via 3d 9 + hv — > 
2p 3d 10 (screening) + e* — >• 3tf 8 + e* + k, where the 3d 10 in 
the intermediate state is due to the screening of the core 
hole. In principle, a Ni:3rf 10 in the ground state would 
be also possible. But we can exclude a strong Ni:3d 10 
based on the XAS data, where we and others [3] observe 
a strong XAS peak at the Ni L3 edge, which would be 
absent for a 3d 10 configuration of Ni. 

The LVV spectrum for Cu is even more structured, 
as demonstrated in Figs.[2]c) and f). Again, we compare 
the experimental results to an atomic multiplet calcula- 
tion [5T] and find excellent agreement with a Cu:3<i 8 final 
state multiplet. This identification is further supported 
by comparing our data to previously published results 



on Cu-metal and CuC>2, which also show the Cu:3rf 8 fi- 
nal state multiplet [24j [25] . 

The Cu AES data therefore unambiguously identify a 
significant Cu:3<i 10 or Cu:3<i 9 component in the ground 
state. In fact, since XAS shown in Fig.[2]c) agrees very 
well with what is expected for a Cu:3cZ 10 [55], we can ex- 
clude a significant population of Cu:3d 9 . This result is 
also in accord with the high (e§ u ) of the Cu:3d states de- 
termined by the PES data in Fig.[l] A Cu:3d 10 state was 
also deduced in a recent x-ray spectroscopy study [1] and 
indirectly concluded from macroscopic measurements on 
BaCu2As2 '-'7 . We can therefore safely conclude that 
Cu is predominantly in a Cu:3rf 10 configuration, i.e., the 
Cu-impurities have a closed shell. 

In the current situation, a realistic estimate can be 
made for the effective onsite Coulomb interaction U e f / on 
the impurities, using the Cini-Sawatzky (CS) theory |28F 
[30] . Within this approach the AES intensity is 



(l-U eff F(e)f + (nU eff D(e)r 

where e is the energy of the 2-hole Auger final state rela- 
tive to the chemical potential, D(e) = J p(x) p(e — x)dx, 
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FIG. 3. Simulated AES spectra and estimated U e ff. Non- 
interacting case: The continuum described by D(e) shifts to 
higher energies with increasing Z due to the increased binding 
energy of the impurity states. Interacting case: while no 2- 
hole bound state is formed for Co, the U B ff is large enough 
to push a 2-hole bond state out of the D(e)-continuum. 



F(e) = P J D(x)/(e — x)dx and p is the impurity DOS 
with the occupied part normalized to 1. Here, we de- 
termined p by means of density functional theory. These 
band structure calculations were done in the local density 
approximation using the full potential code WIEN2K [31], 
as described in Ref.[13l The results of this calculation for 
non-interacting (U e ff — 0, /(e) = D(e)) and interacting 
{U e ff ^ 0) holes are presented in Figj3] 

For Fe (not shown) and Co we determined U e ff by 
fitting the center of mass of the calculated spectra to 
the corresponding experimental values, which are (4.0 ± 
0.5) eV for Fe and (4.5 ± 0.5) eV for Co. For Ni and Cu, 
we set U e f f as to match the measured binding energy of 
the l G term, which is 7.35 eV and 14.85 eV, respectively. 

The energies of the Auger final states are ex- 

pressed on the binding energy scale and were determined 
using the resonant Auger Raman spectra. In this way we 
obtain U eff = (1.4±0.6)eV, (1.8±0.6)eV, (3.0±0.4)eV 
and (7.5 ± 0.4) eV for Fe, Co, Ni and Cu (inset of Fig.|3j. 
?7 e //-values for Fe and Co were also reported previ- 
ously [HI [531 [3D] and are in fair agreement with our re- 
sults. The analysis explains the development of the quasi- 
atomic LVV: with increasing Ueff the Auger final states 
are pushed out of the continuum of states given by D(e) 
and form localized bound states 28 . The emerging mul- 
tiplet structure with increasing Z in Fig. [2] is therefore a 



direct consequence of the increasing U e ff- 

The U e f f determined here corresponds to the effective 
Coulomb interaction between 2 holes in the valence shell. 
U e f f depends on the spatial structure of the 2-hole wave 
function and, most importantly, the screening of these 
holes by the surrounding charges. Especially the latter 
strongly reduces the bare (atomic) Coulomb repulsion 
of ^25 eV to the observed U e ff- The increase of U e ff 
with Z is consistent with the observed charge accumula- 
tion, i.e., the filling of the TM:3d-shell with increasing Z, 
as this indeed reduces the number of possible screening 
channels fcf. Ref. l3"2"]) . It is important to point out that 
the onsite repulsion parameter U used in various theories 
is model dependent, because different approximations for 
the wave functions are used and the screening is taken 
into account to different degrees. U can therefore differ 
from the fully screened U e f / determined here. For exam- 
ple, U for a multi-band model, which takes into account 
all the relevant screening channels, would be close to the 
atomic value of ~25 eV. 

We now discuss the consequences of the above results. 
Our data confirm the strong overlap and hybridization of 
the Co and Fe states. The Co states therefore have a sig- 
nificant Fe:3d character, i.e., Co adds electrons to the Fe 
bands and the number of charge carriers with Fe:3d char- 
acter increases. But this does not mean that the average 
local electron density (n Fe ) at the Fe-sites increases ac- 
cordingly, because the delocalized charge carriers have an 
increased probability of being around Co. Charge density 
therefore piles up around the Co O Q21 [23] , resulting in 
(riFe)— const., as was indeed observed by local probes [9]- 
[TTj . The scattering potential due to the Co- impurities 
is given by (eg°) ~ 1.75 (eg e ) ((eg e ) = 0.4eV)[33] and 

JjCo ^ i o TjFe 
U eff - L - aU eff 

The Ni states are located close to the bottom of the 
Fe bands and the data in Fig. [I] show that these impu- 
rity states contribute only little to the states close to ep. 
At the same time, the Ni:3ci 9 observed here is direct and 
clear experimental proof for the charge accumulation at 
Ni. As in the Co-case, this is perfectly in line with the 
(np e )— const, observed by local probes [9] . The scatter- 
ing potential due to the Ni-impurities is much more se- 
vere than in the Co-case, because (e^ 1 ) ~ 3.5 (e^ e ) and 



U eff 



2.1 Ueff- This will affect the electronic structure 



and result in a broadened FS. 

The Cu impurities differ from the previous cases in 
that they form bound states well below the Fe bands 
with little Fe:3d character. The observed Cu:3<i 10 con- 
figuration implies a strong charge accumulation at Cu, 
which may even result in hole doping of the Fe bands 
(see also Ref. [1]). So the electron doping is reduced and 
possibly even turned into hole doping when going from 
Co over Ni to Cu. We also find that the scattering due to 
Cu with its closed shell is very strong {(e% u ) — 10 (efj e ), 
— 5.3 U^ff)- Correspondingly, a few percent Cu 
substitution will already make the definition of a FS dif- 
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ficult. We note that the above is in line with a very recent 
neutron scattering study [33] . 

To conclude, if and how many electrons the TM sub- 
stitution adds to the Fe bands depends crucially on the 
impurity-host hybridization and the energy difference be- 
tween the host and impurity states. These parameters 
depend on the chosen TM and their effects go beyond 
a mere charge carrier doping. Especially the additional 
scattering potentials, which are quantified here, affect the 
low-energy electronic structure, play a role for the phase 
diagrams and need to be considered in realistic theoreti- 
cal models. Interestingly, the HTS upon Ni substitution 
proves that the HTS in the pnictides is robust against 
strong impurity scattering. 
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SUPPLEMENTARY MATERIAL: 
GROWTH AND CHARACTERIZATION OF 

BA(FE,TM) 2 AS 2 (TM=CO, NI, CU) 




FIG. 4. As grown Ba(Feo.95Nio.os)2As2 single crystals. 



All the crystals used for the present study were grown 
by the self flux technique as described in Refs.fMland [15], 
using FeAs as flux. In case of the Ba(Fei_. c Co :c ) 2 As 2 sys- 
tem, details about the crystal growth, the characteriza- 
tion and the physical properties were published in Ref. 1151 
Here we provide additional details regarding the growth 
and the characterization of the Ni and Cu substituted 
BaFe 2 As 2 single crystals. The pre-reacted precursor ma- 
terials FeAs, Fe 2 As, BaAs and metallic Ni or Cu (TM) 
were mixed, leading to a Ba(Fei_ x TM x ) 3 .iAs3.i compo- 
sition. This composition was used to achieve a homoge- 
neous melt at T — 1463 K. The melt was cooled slowly 
under a temperature gradient in a double-wall crucible 
assembly to obtain large and flux-free single crystals of 
Ba(Fei_ x Ni x ) 2 As2 and Ba(Fei_ :c Cu a ;) 2 As 2 . Two exam- 
ples for the obtained single crystals with typical dimen- 
sions are shown in the Fig|4j 

The high quality of the grown single crystals was as- 
sured by several complementary techniques. From each 
batch, several samples were examined with a Scanning 
Electron Microscope (SEM Philips XL 30) equipped 
with an electron microprobe analyzer for the semi- 
quantitative elemental analysis in the energy dispersive 
X-ray (EDX) mode. Using EDX, the Ni and Cu concen- 
tration was determined by averaging over several differ- 
ent points on the sample surface. The estimated com- 
position from the EDX for the present single crystals is 
Ba(Feo.9 5 Nio.o 5 )2As 2 and Ba(Fe .9iCu .o9)2As 2 . In gen- 
eral, the error of an EDX analysis is about 2 mass per- 
cent without any additional standardization procedure. 
However, the method of averaging reduces the size of the 
error bars significantly. In this way the concentration x 
of Co, Ni and Cu was determined to be (0.06 ± 0.01), 
(0.05 ± 0.01) and (0.09 ± 0.01). 

In addition to this, the TM contents of the studied 
samples were also verified in-situ by XPS measurements, 
which gave the same results as EDX. The in-situ XPS 
measurements also do not show any sign of a signifi- 
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FIG. 5. XPS survey scans. Only a small trace of oxygen is 
observed, which can at least partially be attributed to adsor- 
bates on the sample surface and verifies the high purity of our 
samples. 



cant oxygen contamination in any of the studied sam- 
ples. Representative XPS survey scans are displayed in 
Fig. [5] The weak O Is line is comparable to what has 
been reported previously in the literature for high-purity 
samples [1] , which verifies the high quality of the studied 
single crystals. Our XAS measurements also exclude the 
presence of Fe-oxide impurities, which would leave a clear 
fingerprint in these spectra |10j . 

Prior the synchrotron experiments, the high qual- 
ity of our single crystals was also checked by XRD, 
rcsisitivity and mangnetization measurements. In 
Fig. [6] we show representative measurements for our 
Ba(Fe .95Ni .o5)2As 2 single crystals. The data in Fig.[6j 
reveals that the structural and magnetic transitions are 
completely suppressed by 5% Ni-substitution, whereas 
bulk superconductivity occurs. This is in excellent agree- 




T(K) 

FIG. 6. Magnetization and resistivity measurement of 
Ba(Fe .95Nio.o5)2As 2 single crystals 
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ment with the literature (S.L. Bud'ko et at, Phys. Rev. 
B 79, 220516R (2009)). Furthermore, the superconduct- 
ing transition at Tc = 21 K is very sharp and, within 
the error bars of our experiment, we find a supercon- 
ducting volume fraction of 100%. This further confirms 



the very high quality of our single crystals. Also the 
Cu-substituted samples were characterized by the same 
methods prior the synchrotron experiments. The struc- 
tural quality was verified by XRD and, in agreement with 
the previous literature, the sample with 9% Cu substitu- 
tion does not show superconductivity down to 1.8 K. 



